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In order to explore how superconductivity arises when charge fluctuations and spin fluctuations 
coexist, we have obtained a phase diagram against the off-site repulsion V and band filling n for the 
extended, repulsive Hubbard model on the square lattice with the fluctuation exchange approxima- 
tion. We have found the existence of (i) a transition between dxy and dx2_y2 pairing symmetries, (ii) 
/-pairing in between the dx2_y2 and CDW phases for intermediate 0.5 < n < 1.0 and large V, and 
(iii) for anisotropic cases the pairing symmetry changing, in agreement with a previously proposed 
"generic phase diagram", as d ^ f ^ s when V (hence the charge fluctuations) are increased. All 
these are consequences of the structure in the charge and spin susceptibilities, which have peaks 
habitating at common or segregated positions in k space. 

PACS numbers: 74.20.Mn 



I. INTRODUCTION 

It is gradually becoming clearer that, while spin fluctu- 
ations are usually considered for electron mechanisms of 
high-Tc superconductivity, charge fluctuations may pos- 
sibly play essential roles in some of the unconventional 
superconductors. Among these are organic metals such 
as (9-(BEDT-TTF)2Xi or (BEDT-TTF)3Cl22H202 that 
exhibit superconductivity sitting adjacent to the charge 
density wave (CDW) in the phase diagram. This suggests 
that charge fluctuations can mediate pairing as well, just 
as the spin-fluctuation-mediated pairing^ can appear ad- 
jacent to the spin density wave(SDW) phase as in the 
cuprates and some of the organic superconductors like 
/^-(BEDT-TTF)2X. 

Studies for the pairing mediated by charge fluctuations 
have been rather scant. Scalapino et al. have shown, 
with the random phase approximation, that dj.2_y2 pair- 
ing gives way to dxy in a three-dimensional cubic lattice 
when charge fluctuations become large with the intro- 
duction of the off-site interaction.^ However, systematic 
studies are yet to come for the charge fluctuations, in con- 
trast to the spin-fluctuation mediated superconductivity 
for which favorable situations for its occurrence has been 
extensively discussed^ 

There is in fact one proposal in the context of the spin- 
triplet superconductivity in a quasi-one-dimensional (ID) 
organic metal (TMTSF)2PF6: Three of the present au- 
thors have proposed a "generic phase diagram" in which 
the dominant pairing symmetry changes sls d ^ f ^ s 
as the charge fluctuation becomes stronger^. The physi- 
cal background is as follows: triplet pairing is very hard 
to be realized to start with, for the reasons identified 
in Ref.^, where the primary one is the strength of the 
pairing interaction for triplets being only 1/3 of that for 
singlets. The situation can be reverted when charge fluc- 
tuations are dominant as discussed in Kuroki et alJ ^ but 



the charge fluctuation was treated phenomenologically 
there, so a microscopic study is highly desirable. Namely, 
while spin fluctuations dominate over charge fluctuations 
when the electron-electron interaction is short-ranged (as 
in the Hubbard model), charge fluctuations should be- 
come more intense as we increase the range of the inter- 
action. The problem becomes especially intriguing when 
charge and spin fluctuations coexist, since they may in- 
duce quantum phase transitions among different pairing 
symmetries. 

This is exactly our motivation here to study the ef- 
fect of strong charge fluctuations by adopting the ex- 
tended Hubbard model, as a simplest one in which we 
can control the relative magnitude of the charge fluctu- 
ation by varying the off-site Coulomb repulsion V. The 
extended Hubbard model has been studied, primarily for 
specific charge densities n, e.g., half- filling or quarter fill- 
ing, by means of quantum Monte Carlo method^, weak 
coupling theory^, mean- field approximation^^ second- 
order perturbation^^, random phase approximation^, 
fluctuation exchange (FLEX) approximation^^, slave- 
boson technique^^, bosonization and renormalization 
group^^'^^. However, the phase diagram of the two- 
dimensional extended Hubbard model against n and V 
has yet to be obtained. 

Here we have determined the symmetry of the dom- 
inant pairing in the V-n space for the extended Hub- 
bard model by focusing on the case of isotropic or 
anisotropic square lattice, since many of unconventional 
superconductors are two-dimensional(2D) or quasi-one- 
dimensional(lD). We adopt the FLEX developed by 
Bickers et al}^~'^^ ^ which is a renormalized perturbation 
scheme to study pairing instabilities when exchange of 
spin and charge fluctuations are considered as dominant 
diagrams. Although this is an approximation, we can 
explore the tendencies for pairing when the system pa- 
rameters are varied. 

We find that (i) there exists a phase transition be- 



tween dxy and dx2_y2 pairing symmetries, (ii) triplet /- 
pairing appears in between the dx2_y2 and CDW phases 
for intermediate 0.5 < n < 1.0 and large y, and (iii) for 
anisotropic cases the pairing symmetry changes, in agree- 
ment with the generic phase diagram^, as d ^ f ^ s 
when V (hence the charge fluctuations) are increased. 
Origin of all these has been identified as the structure 
in the charge and spin susceptibilities, which can have 
peaks habit ating at segregated positions in k space. 



II. FORMULATION 

Let us start with the extended Hubbard Hamiltonian, 

nn 

W = - ^ ^ Ujcl^Cj^ + ^ XI ^^T^a 

i,j a i 
^ nn 



in the standard notation on a tetragonal lattice de- 
picted in Fig. 1. For the (isotropic) square lattice 
the unit of energy is taken to be the nearest-neighbor 
tij = 1.0, and lattice constant a = 1. 
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FIG. 1. A tetragonal lattice with hopping integral tx 
along X-axis and ty along along with the nearest-neighbor 
Coulomb repulsion Vx along x-axis and Vy along y. 



To determine the dominant gap function, we solve 
Eliashberg's equation with the FLEX approximation. 



(2) 



where is the gap function, G Green's function, and F 
the pairing interaction with k = {k^uOn)- The eigenvalue 
A, a measure of the pairing, becomes unity at T = Tq. 
For the calculation we take an = 32 x 32 lattice, the 
temperature T = 0.02, and the Matsubara frequency 
for fermions -(2 Ac - 1)t:T < uJu < {'^.Nc - 1)t^T with 
Ac = 1024. 

Esirgen et a/.^^'^^'^^ have extended the FLEX method 
to general lattice Hamiltonians including the extended 
Hubbard model. Following them we introduce the pair- 
ing interaction. 



Ts{k,k') = 

J2 { i[KnXspK.](fc-fc';Ar;Ar')e^(fc-^^+^'-^^') 



Ar,Ar' 



+ i V;{0; Ar; A^')e'(fe-Ar'-fc'.Ar) } ^ (3) 
for singlet pairing, and 

Tt{k,k') = 

E {- H^mXspK.](fc-fc';Ar;ArV('^-^^+^'-^^'^ 



Ar,Ar' 



- \ [V^Xci.V^]{k-k';Ar;Ar')e'(k-Ar+k'-Ar') 
+ \ yt(0; Ar; ^r')e'^k-Ar' -k' -Av) | (4) 



for triplet pairing. Here Ar(= 0,=bx,it^) is null or 
nearest-neighbor vectors, 

Xsp =X/(l + KnX), 

Xch = x/(l + KiX) 

are the spin and charge susceptibilities, respectively, 
where x is the irreducible susceptibility, 

xfeAr;Ar') = ^ e'^'-(^^-^^')G(fc' + 5)G(fc'), (5) 



and Vd(14i) is the coupling between density (magnetic) 
fluctuations, 

yd(g; Ar; Ar ) 

r U + 4[14 cos(^,,) + Vy cos(^^)], Ar = Ar' = 0, 
= \ -Vx, Ar = Ar' = ±x (6) 

[ -Vy, Ar = Ar' = ±y 

( -U, Ar = Ar' = 0, 
V^{q- Ar; Ar') = <^ -K, Ar = Ar' = ±x (7) 
[ -Vy, Ar = Ar' = ±y, 

where q = {q^en) with = 2n7rT being the Matsubara 
frequencies for bosons. We have found that the q de- 
pendence of Vjn and Vd does not in fact affect Fg and Ft 
significantly. Accordingly the peak position of Xch is al- 
most the same as that for VaXchKi term in the expression 
for F. 

14(0; Ar; Ar'), Vt{^] /^r] /^r'), appearing in the last 
lines in eqs.(3,4) respectively, are constant terms involv- 
ing U and y, 

2U, Ar = Ar' = 0, 
14, Ar = Ar' = ±x 
Vs{q; Ar; Ar') = { y,,e±^^-, Ar = -Ar' = ±x (8) 

Vy, Ar = Ar' = ±y 
Vye^'^y, Ar = -Ar' = ±y 

Vx, Ar = Ar' = ±x 
A A /X ) -Vxe^'"^-, Ar = -Ar' = ±x 
y,(g;Ar;ArO=<! yj^ Ar = Ar' = ±y 

-VyC^'^iy, Ar = -Ar' = ±y 



When the oflF-site interaction V is introduced ah the 
vertices (Fm, Vd, Vs, Vt) as weh as the susceptibihties be- 
come (Z+1) X (Z+1) matrices for the lattice coordination 
number Z(= 4 for the square lattice). 



III. RESULT 



A. Square lattice 

Let us first display the obtained phase diagram against 
V and the band filling n for the square lattice in Fig. 2. 
The phase diagram is drawn by assuming that the pair- 
ing instability that has the largest A in Eliashberg's equa- 
tion, calculated at T = 0.02 here, has the highest transi- 
tion temperature, because, with the present complexity 
of the model, it is difficult to extend the FLEX calcu- 
lation to lower temperatures. While the value of A for 
T = 0.02 is still much smaller than unity (see below), it 
is expected that the order in which the values of A for 
various phases appear do not change for T ^ 0. The 
density waves (CDW or SDW) are identified as the re- 
gion where the respective (charge or spin) susceptibility 
diverges (at T > 0.02 in the present calculation). The 
area of the CDW and SDW phases may expand at lower 
temperatures, but the phase diagram is not expected to 
change significantly. 

Figure 2 is thus obtained, with the on-site Coulomb 
repulsion fixed at = 4 hereafter, and we immediately 
note that CDW, SDW, singlet superconductivity (SC), 
and triplet SC all appear on the (V, n) plane. The phase 
diagram is reminiscent of that for one-dimensional ex- 
tended Hubbard model obtained with the Tomonaga- 
Luttinger theory However, an essential difference is 
that U or V has to be negative (attractive) to realize SC 
phases in one dimension, while we are talking about the 
case when both U and V are repulsive in 2D. Another 
comment is that in our calculation we cannot treat the 
Mott insulator, which should appear at n close enough 
to the half-filling (n = 1). 

Before elaborating on the superconducting phases, let 
us make a remark on density waves. Intuitively, CDW 
should appear for strong enough V, while we should have 
SDW for the band filling close enough to the half-filling. 
The boundary between CDW and SDW for n ^ 1 is 
seen in the present result to fall upon a line representing 
y = 1, which agrees with a mean-field picture: In the 
CDW state where electrons doubly-occupy every other 
sites, each electron feels on average an on-site energy 
U/2 per electron, while in the SDW state each electron, 
singly-occupied, feels an off-site repulsion V/2x Z = 2F, 
so the SDW/CDW boundary corresponds to V = 
As for the wave vectors describing the CDW (SDW) or- 
dering, they are as indicated by the peaks (around (tt, tt)) 
in the charge (spin) susceptibilities, which are displayed 
below. 
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FIG. 2. Phase diagram against V and n with U/t = 4 for 
the 2D extended Hubbard model. 



A salient feature for superconducting phases in Fig. 
2 is that a spin-triplet /-phase appears just below the 
CDW phase for an intermediate region of n and V. The 
behavior of A when V is varied with a fixed n = 0.7 is 
depicted in Fig. 3.^^ Figure 4 shows the gap function in k 
space for the /-wave, which should be called, more pre- 
cisely, in the group theoretical representation. This 
has only one nodal line passing F point, but may be 
called / in that the gap function (~ sm{kx) + sm{ky) + 

const. [sin ( 2 /ca;) +sin(2A:^)]) changes sign as H \ \ — 

around the Fermi surface. 




FIG. 3. The maximum eigenvalue, A, of Eliashberg's equa- 
tion for the spin-singlet (solid line) and triplet (dotted) chan- 
nel as a function of 1/ for n = 0.7 with the dominant orbital 
symmetry indicated. The CDW phase is identified from the 
divergence in the charge susceptibility. 
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FIG. 4. The /-wave (irreducible representation F^) gap 
function in k space (right panel) for n = 0.7 and V = 1.3, 
as compared with the (1^2 _y2 (left) for n = 0.7, V = 0.5. 
Nodal lines are represented by dotted lines, while the Fermi 
surface by solid curve on which the sign of the gap function 
is indicated. The arrow indicates a typical scattering process 
mediated by spin fluctuations. 



FIG. 5. Top: the maximum eigenvalue of the charge (solid 
curve) and spin (dotted) susceptibilities as a function of V for 
n — 0.7. Bottom: color-coded plots of the charge (left) and 
spin (right) susceptibilities in the (kx,ky) space for V = 1.3 
and n = 0.7. The susceptibilities (energy)"^] in units in 
which t = 1. 



B. Physical origin — common vs segregated peaks 

in Xsp and Xch 



We can keep track of the origin of the spin-triplet 
instability by looking at the structure (peak intensi- 
ties and peak positions in k space) of the charge (Xch) 
and spin (Xsp) susceptibilities (static with Matsubara 
frequency = 0) in Fig. 5. The peak intensity of Xch is seen 
to exceed that of Xsp as V is increased. If we turn to the 
structure in k space, Xch and Xsp have similar peak po- 
sitions. If we go back to eqs.(3),(4), the spin- fluctuation 
mediated pairing interaction (the first line on the right- 
hand side) and the charge-fluctuation mediated one (the 
second) act destructively (with opposite signs for the two 
terms in eq.(3)) for singlet pairing, whereas they act con- 
structively for triplet pairing (same signs in eq.(4)). So we 
can expect the realization of triplet pairing (|rt| > iFgl) 
for large enough charge-fluctuation mediated pairing in- 
teractions (i.e., for large enough V) provided Xch ctnd Xsp 
have common peak positions. 

Physically, Xsp = x/(l - ^x) while Xch = x/(l + ^x) 
when y = 0, so the nesting of the Fermi surface (which 
dominates x) determines the peaks of Xsp (but not those 
of Xch)- When V is switched on, Xsp = x/(l + ^mX) and 
Xch = x/(l + ^x) ^ay have common peak positions, 
which is what is happening here. 




How the situation is altered for smaller n for which dxy 
pairing turns out to be dominant? In Fig. 6, which shows 
Xch and Xsp for ^ = 0.6 with V = 1.6, Xch has a larger 
peak than Xsp does, but the peak positions of Xch are dis- 
tinct from those of Xsp- Namely, while Xch has peaks at 
four points around Q^j^ :^ (±7r, ±7r), Xsp has peaks that 
are shifted toward Q^^ = (±7r, 0) and (0, ±7r). Now, as a 

basic property of Eliashberg's eq.(2) with eq.(3) plugged, 
to realize a large A requires that (i) the gap function (j) 
must change sign (across the typical pair-scattering mo- 
mentum transfer Qgp) to turn the originally repulsive 
spin- fluctuation- mediated interaction (i.e., the VmXspVm 
term with a positive coefficient, 3/2 in eq.(3)) into an 
effective attraction in the gap equation, while (ii) <p must 
not change sign (across the pair-scattering momentum 
transfer Q^h) make the originally attractive charge- 
fluctuation- mediated interaction (the FdXch^d term with 
a negative coefficient, —1/2) remain attractive. 

We can see in Fig. 7 that dxy pairing does satisfy the 
above condition. For n closer to half- filling, by contrast, 
Xsp becomes dominant with Q^^ (i^r, ±7r) (an arrow 
in Fig. 4), so dx2-y2 is favored as usual. So the picture ob- 
tained here is that the pairing symmetry can change, even 
within the spin-singlet channel, when common peaks be- 
tween Xsp and Xch change into segregated peaks (as n 
and/or V are changed). This is contrasted with the above 
case of triplet /, for which a transition between singlet 
and triplet {d ^ f) can occur with the common peaks 
between Xsp and Xch throughout. For n = 0.5, Merino et 
a/.^^ and Kobayashi et al.'^^ have recently derived similar 
results using the slave-boson technique and RPA, respec- 
tively which are consistent with the present result. 
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FIG. 6. A plot similar to Fig. 5, for n = 0.6 with V = 1.6. 
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FIG. 7. The Fermi surface for n — 0.6, V = l.Q with dxy 
pairing (right), and for n = 0.7, V = 1.3 with the / pair- 
ing(left). The black (gray) arrows indicate typical scattering 
processes mediated by charge (spin) fluctuations. 



C. Quasi one-dimensional lattice 
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FIG. 8. Result for a quasi ID system with tx = 1.0, ty = 0.2 
for n = 0.9. Top: a plot similar to Fig. 3. Bottom: The dom- 
inant gap function (with nodal Unes represented by dotted 
lines) for K = 1.0 {d wave), Vx = 2.2 (/), and Vx = 2.4 (s), 
along with the Fermi surface (solid curves) on which an arrow 
represents the nesting vector. 



IV. CONCLUSION 



Let us finally discuss the anisotropic (quasi ID) case. 
While the original proposal made in Ref. for the triplet 
'/-wave' pairing in (TMTSF)2X (where the symmetry 
refers to the warped Fermi surface; see Fig. 8 below) is 
for a quarter- filled case, the competition between the 
charge and spin fluctuations should become more strin- 
gent near half-filling, where U introduces 2kF spin fiuc- 
t nations while V enhances 2k f charge fluctuations. Thus 
the triplet pairing should dominate over singlet d beyond 
some value of V. Note that the triplet gap function has 
to be / rather than p, since, with Tt{Q) negative (at- 
tractive), the gap has to have the same sign across the 
nesting vector Q. When V becomes even larger, rs(Q) 
turns negative (attractive), so that the singlet s with no 
nodes on the Fermi surface will take over. 

In order to show that this prediction is indeed real- 
ized, we have performed a FLEX calculation on a quasi 
ID lattice with tx = 1.0, ty = 0.2. To represent the quasi 
ID system we have here assumed that the off-site re- 
pulsion only acts between nearest neighbors along x (the 
conductive direction), although this assumption does not 
qualitatively alter the result. Figure 8 plots A as a func- 
tion of 14 for n = 0.9, along with the forms of (i, /, and 
s gap functions. We can see that the dominant pairing 
changes 8iS d ^ f ^ s with V, in an exact agreement 
with Ref.^. 



We have obtained the phase diagram for the 2D ex- 
tended Hubbard model with the FLEX approximation. 
We have found that /-wave pairing is favored for interme- 
diate n and large V. In the more dilute regime (n ~ 0.5) 
the peak position of Xch and Xsp is separated, and dx2-y2 
pairing gives way to dxy These are the consequences of 
how each pairing symmetry can exploit charge and spin 
fluctuations, and we end up with a picture that the pair- 
ing symmetry can change, even within the spin-singlet 
channel {dx2-y2 ^ dxy)^ when common peaks between 
Xsp and Xch change into segregated peaks (as n and/or 
V are changed), which is contrasted with the case of com- 
mon peaks throughout with a transition triplet f ^ d. 
In the anisotropic (quasi-one-dimensional) case, the dom- 
inant pairing gap function changes as d f s with 
y, which agrees with the phenomenological theory^. 

In a broad context it should be interesting to examine 
how V (hence the charge fluctuation mediated interac- 
tion) can dominate the pairing symmetry in unconven- 
tional superconductors, although the nature of the pair- 
ing may be sensitively affected by the underlying band 
structure of each material. 
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